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SMALL SCALE TURBULENCE IN THE CRAB NEBULA; 
EVIDENCE OF LOWER HYBRID PARAMETRIC INSTABILITIES 
DRIVEN BY THE PULSAR WAVE 

ABSTRACT 

Strong small scale turbulence with size a ~ 50 km is observed in 
the Crab Nebula from the temporal pulse broadening data. It is 
shown that the strong 30 Hz pulsar wave can parametrically 
excite instabilities near the lower hybrid frequency in the thermal 
plasma of the Crab Nebula with a characteristic wavelength of the 
order of the scale size a of the turbulence observed. The growth 
rate of the excited waves is about 0. 2 - 10 sec. These instabilities 
provide a coupling mechanism between the pulsar wave and the 


Nebula plasma. 


SMALL SCALE TURBULENCE IN THE CRAB NEBULA: 


EVIDENCE OF LOWER HYBRID PARAMETRIC INSTABIUTIES 
DRIVEN BY THE PULSAR WAVE 


(I) INTRODUCTION 

Recent observations of angular broadening and temporal broadening of the 
Crab Nebula pulsar’s radiation indicate the existence of a scattering region of 
plasma turbulence within the Crab Nebula (Counselman & Rankin 1972, Rankin 
& Counselman 1973, Lyne & Thorne 1975, Vandenberg 1976). Ui^mg the results 
of VLBI (very-long-baseline interferometer) data, Vandenberg (1976) was able 
to show that the angular broadening and the temporal broadening of the pulsar 
radiation occur in different locations along the line of sight to the Crab pulsar. 
Vandenberg (1976) then proposed the electron density fluctuation in Filament 
159 of the Crab Nebula (Trimble 1970) as the cause of the temporal pulse 
broadening. However, Vandenberg (1976) made a mistake in estimating the 
scale size and the amount of plasma irregularity in the filament. Lee (1977) 
developed a thin screen scintillation theory for a spherical wave and corrected 
Vandenberg's errors. Lee (1977) pointed out that if the pulse broadening is due 
to the thermal plasma in the Crab Nebula (of thickness ■>'0. 5 pc), then the rms 
electron density fluctuation AN in the Nebula is 

AN/^>„s=0.01 cm'^ (1) 

and the scale size of electron density fluctuation, a, is given by 
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■/^^«lSkm (2) 

where 0 ^ is the mts phase fluctuation caused by the Neintla. From the pulse 
shape observed, one must have In the next section, we estimate 0 ^ to be t 

~l-5. Thus AN ~ 0.01-0.05 cm~^ and a ~ 15-75 km. However, if the pulse 
broadening is due to Filament 159 (with thickness ~0.01 pc) as proposed by 
Vandenberg, then one has in the filament AN/0„ 0. 08 cm“^ and a/0^j ~ 15 km. 

It is the purp<»se of this paper to propose a mechanism to explain the ob- 
served small scale turbulence in the Crab Nebula. We will find in the next 
section that the averaged thermal electron density in the Nebula is N^^^ 0.05 

cm“^ . Thus, for the small scale turbulence (a «« 50 km), AN ~ Nj,j and the 
turbulence is very strong. We believe that some external driving force is 
needed to excite such a strong turbulence at a relatively small scale size. The 
natural candidate for the driving force is the strong 30 Hz pulsar wave. In 
Section 11, we briefly describe the physical environment in the Crab Nebula. In 
Section III, we show that the 30 Hz pulsar wave can indeed excite several para- 
metric instabilities near the lower hybrid frequency in the Crab thermal plasma 
with a characteristic wavelength a. Thus these instr.bilities can explain the 
observed strong turbulence and, in addition, can provide a coupling mechanism 
between the pulsar and the nebula. In Section IV, we discuss and sum up the 


results. 
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ai) NEBULA 

(a) THERMAL PLASMA 

The Crab Nebula contains three components; relativistic charged parti- 
cles, magnetic field and thermal plasma. The magnetic field in the Crab 
Nebula has been observed to be ~ 3 X 10“^ Gauss (e. g. , Fazio et al. 

1972, Fazio 1973) and the total energy of relativistic electrons is about 10^^ 
erg (Trimble & Rees 1970). The thermal electrons in most filaments have a 
density, N,^ ~ 10^ cm*^ and a temperature, T^ 10^ °K (Woltjer 1958, 
Trimble 1970,. The thermal plasma density outside filaments is much lower 
than that in filaments. Shklovskii (1966) pointed out that in the Crab Nebula 
the averaged thermal plasma density, N^^ < 1 cm"^ , because of the absence 
of the gross Faraday depolarization of radio emission. From the motion of 
wisps, Scargle (1969) argued that N^j^ ~ 5 x 10"^ cm“^ near Wisp 1 and 
Nj^ ~ 0. 05 cm'^ at the northwest edge. From the fact that the expansion 
speed of the nebula at the outer edge, ~ 1,5 X lO* cm /sec, much less 
than the speed of light, one would obtain N^^ = 0.05 cm"^ for * 10’* 
Gauss near the outer edge of the nebula (Kahn 1971). Assume that N,^ > AN 
(Lee & Jokipii 1976). From Eq. (1), we also find that N,^ > 0.01 cm'^. Thus 
we take N^j^ ~ 0.05 cm'^ . If AN < N^jj as assumed, Eq. (1) gives 1-5. 
Therefore, AN ~ 0. 01-0. 05 cm"^ and a ~ 15-75 km. Assuming energy equi- 
patition between magnetic field and thermal plasma in the Nebula, we have the 
thermal plasma temperature T^ ~ 5 x 10^ ®K, Thus in the following 
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calculations, we will take =r 3 X lO'^Gauss, N,jj= 0.06 cm"^ and 
T ~ 5 X 10'' "K. 

e 

(b) 30 Hz PULSAR WAVE 

Crab pulsar is generally believed to be a compact oblique magnetic 
rotator and therefore emits strong 30 Hz electromagnetic waves (Ostriker & 
Gunn, 1969). Just outside the speed-of-light cylinder, the wave electric field 
E is strong and the parameter 


where c is the speed of light, e the electron charge, the electron mass 
and = 30 Hz is the frequency of the pulsar wave. This 30 Hz wave can 
accelerate particles, leaking from the pulsar magnetosphere, to relativistic 
velocity (Ostriker & Gunn 1969, Ferrari & Trussoni 1975). The lack of cir- 
cular polarization in the optical continuiun radiation from the nebula strongly 
suggests that the pulsar wave does not propagate freely more than ~20" (or 
6 X 10*’ cm) from the pulsar (Landstreet & Angel 1971, Rees 1971). Observa- 
tions of a 'valley' in the brightness near the pulsar at both radio and optical 
frequencies suggest there exists a cavity around the pulsar swept clear of 
thermal gas and static magnetic field by the 30 Hz wave (Scargle 1971). 

Thus we have a picture for the Crab Nebula: the 30 Hz pulsar wave 
propagates from the pulsar and interacts, at a distance of several 10*’ cm 
away from the pulsar, with the Crab Nebula, which is composed of t'nermal 
plasma and magnetic field. At the edges of the central pulsar cavity, the wave 
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pressure becomes comparable to the thermal and magnetic pressure of the 
Nebula and the wave field E « . As the wave propagates into the Nebula, 

the pulsar wave will be damped and/or reflected such that E < B ^inside the 
nebula. 

Several authors have considered the Interaction between the pulsar wave 
and the nebula. Rees & Gunn (1974) suggested that the 30 Hz wave would be 
absorbed due to synchrotron absorption by the downstream shocked high energy 
plasma with a mean free path ~3 X lo'^ cm. Ferrari (1974) and Dobrowolny 
& Ferrari (1976) discussed the plasma instabilities excited by the interaction 
between the 30 Hz pulsar wave and the relativistic wind leaking from the pulsar 
magnetosphere. These instabilities occur only near the equatorial plane of the 
pulsar spin axis where wisps are observed. Max (1973) discussed parametric 
instabilities excited b> the 30 Hz wave in the nebula thermal plasma, neglecting 
effects of the nebula magnetic field. The existence of a static magnetic field 
allows more possible modes to be excited. In the next section, we will consider 
the interaction between the 30 Hz wave and the thermal plasma in the presence 
of the nebula magnetic field. 

(Ill) PARAMETRIC EXCITATION BY THE PULSAR WAVE 

In the Crab Nebula where N,jj ~ 0. 05 cm*^, B^^j, ~ 3-X 10'^ Gauss, and 
Tp ~ 5 X lO’ °K, we have the electron cyclotron frequency 840 Hz, ion 
cyclotron frequency I2. ^ 0.5 Hz, electron plasma frequency 1800 Hz, 

ion plasma frequency ~ 40 Hz, Debye length Xp ~ 2 km, and the ion sound 



speed V ~ 8 X 10^ cm/sec. The lower hybrid frequency w, „ is given by 

‘^LH = ‘^pi ^ 20 Hz (3) 

Thus the 30 Hz pulsar wave, when propagating in the thermal plasma of the 
Crab Nebula, acts as a pump wave with frequency in the vicinity of the 
lower hybrid frequency ( ^ Wj_jj ). Parametric decay instabilities with 
pump frequency near has recently been studied (Kindel et al, 1972, Porkolab 
1974, Hesegawa & Chen 1975, Better b Perkins 1976, Ott 1975, Tripathi et al. 
1977) and observed in the rf heating of tokamak plasma (Hooke b Bemabei 1972, 
Chu et al, 1973, Chang b Porkolab 1974, Porkolab et al. 1977). 

Consider now the i^irametric decay of the pump wave (with - 30 Hz 
and wavenumber kp) Into two daughter waves ( Wj , jc,) and ( Wj, kj). Since the 
wavelength of the incoming pump wave is much greater than that of the daughter 
waves (~10^ cm), the pump will be treated in the dipole approximation (k^ 0), 

£= cosfco^t). (4) 

We then have kj k.2 “ !So Under the dipole approximation, Porkolab 
(1974) has investigated possible decay channels for a pump wa^'e near lower 
hybrid frequency. According to Kindel et al. (1972) and Porkolab (1974), the 
pump wave will excite (a) purely growing and (b) decay parametric instabilities. 
In both cases, one of the two daughter waves is a lower hybrid wave. Let 
(coj, k,) be the excited lower hybrid wave satisfying the following dispersion 


relation 
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«2=WLk = ‘*'LH 


" kl M j ’ 


(5) 


where w, „ the lower hybrid frequency is given by Eq. (3), M, is the proton 

LH •'I 

mass and kjn is the component of the wavevector along the magnetic field, 
(a) PURELY GROWING PARAMETRIC INSTAEILITIES 


For the purely growing instabilities which occurs for < Wj , the real 
part of cjj is zero and the minimum threshold is given by 


“p.U\ "'/J 


( 6 ) 


where U = CE /B . is the E X B drift velocity, V = (kT /m is the ion 
sound speed and rj is the linear damping rate of the lower hybrid wave. The 
damping rate for the lower hybrid wave is 


2 - ,.,2 . 


=V7- 


- UJ' 


LH 




y^exp(- y^) 


(7) 


where y s w/k||Vj^ and is the electron thermal velocity (Tripathi et al. 
1977). Since y ^ 1, ■)'2 very small and therefore the condition In Eq. (6) is 
easily satisfied for the pulsar wave propagating in the Nebula. The growth rate 
for the purely growing instabilities well above threshold is given by 

7~(V^X 

where |3 - ( ‘*'2' ’^®*ting x = 1, = 10"^ for 

example, we have y-z 0.07 ^1.4 sec'* . In general the t.me scale for 

the purely growing instaHilities in the nebula is of order of 1 sec. 
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(b) DECAY PARAMETRIC IN8TA0IUTIES 

For decay instabilities, the pump wave beating with a lower hybrid wave 
form low-frequency ion quasimodes with frequencies near or above the ion 
cyclotron frequency (Porkolab 19741. These quasimodes may be of two types: 

(i) kinetic modes such that their parallel phase velocities are near the electron 
thermal velocity ( o>/k|| ^ which case the driving mechanism is essen- 

tially nonlinear electron Landan damping; (ii) nonresonant fluid-like modes such 

that uj/k n ^ V . 

Ii f e 

For kinetic modes, the threshold is the same as that for the purely grow- 
ing instabilities in Eq. (6). Let w, = ^ iU| . The growth rate well 

above threshold is 


t0|/w2 i 



Vny exp (-y^) 


(9) 


where y s w^/kiiVj^ ~ 1. However, from the numerical calculation of 
Porkolab (1974), co, is peaked around 1. The maximum growth rate 

can then be obtained from Eq. ( if) by setting U/V^ = 1, cj|(max) :r 0.02 wj 
0. 6 sec ‘ . In general the time scale for the growth of the kinetic modes is 
about 1-10 sec. 

The fluid-like quasimodes ( tj^/k(| > ) can easily be excited even in 

weak pump fields (Porkolab 1974). The growth rate is given by 
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Again the growth rate w, is peaked at U/V, 1. From Eq. (10), the maximum 
growth rate is w,(max) ir 0. 2 =: 6 see’* . Thus the time scale for fluid* 

like modes in the nebula is about 0.2-5 sec. 

The wavelength of the excited waves for both purely growing and decay 
parametric instabilities is given by kX^ - 0. 1-0. 5 (Porkolab 1974, TrifXithi 
et al. 1977). In the Crab Nebula, we have 2 km, and therefore the wave- 
length X = 2)r/k ~ 20-120 km, which is approximately the scale size a (—15-75 
km) of the turbulence observed. 

(IV) DISCUSSIONS 

In the above calculation, we assume that the 30 Hz pulsar wave reaches 
the lower hybrid resonant region directly. If the thermal plasma 

density increases gradually near the inner edge of the Crab Nebula, then before 
reaching the resonant region, the 30 H’' pulsar wave will pass through a layer 
where and may be reflected. However, since the pulsar wave is 

strong (a > 1), it can penetrate through this layer (Max and Perkins 1971) and 
reach regions where - ^Ltr whereas in the tokamak plasma 

the lower hybrid resonant region is only a thin layer, the lower hybrid resonant 
region in ihe Nebula is very broad since in the Nebula is determined 
mostly by the magnetic field ( 1*^ ), not by the plasma density 
(or u)pj). 

Next we consider if the temporal pulse broadening is due to the thermal 
plasma fluctuations in Filament 159 as proposed by V^andenberg tl976). In 
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Filament 159 where 10^ °K and ~ 100 cm"^ , the Debye length is 
Ap ~ 70 cm. If waves are parametrically excited in the filament by the pulsar 
wave, then the wavelength A would be only ~(l-4) X lO^ cm which is too small 
to account for the observed turbulence. 

Thus the .30 Hz pulsar wave can excite parametric instabilities near the 
lower hybrid frequency in the thermal plasma of the Crab Ne'jula w th a charac- 
teristic wavelength of the order of the scale size a (~15-75 km) of the turbulence 
observed. The time scale for the growth of the excited waves is about 0. 2-10 
sec. The excited waves propagate nearly perpendicular to the Nebula field 
(k||< k) and therefore will propagate radially outward with a phase velocity 

10® cm/sec superposed on the convective motion of the Nebula plasma. 
Since the damping rate of the lower hybrid waves as given by Eq. (7) is very 
low, these waves will propagate to the outer part of the Nebula. We note that 
as in the rf toxamak plasma heating, the electron.^ and ions in the nebula will 
ixj heated by this mechanism. If we identify the scale size a of the observed 
turbulence with the wavelength of the excited waves (a ~ A 2.. * ^ (0. 1-0, 5) x 
Ap), then A,j gives the plasma temperature in the Nebula, - 10’ -10* "K, 
which is consistent with that determined by energy equipai iition in Section 11. 
Thus the lower hybrid instabilities also provide a coupling mechanism lx;tween 
the :)0 Hz wave and the surrounding Nebula, in which the energy and momentum 


of the pulsar wave 


transferred to the plasma in the Nebula. 
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